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We present the numerical solution of the renormalization group (RG) equations derived in RefQ, 
for the problem of superconductivity in the presence of both electron-electron and electron-phonon 
■ coupling at zero temperature. We study the instability of a Fermi liquid to a superconductor and the 

RG flow of the couplings in presence of retardation effects and the crossover from weak to strong 
coupling. We show that our numerical results provide an ansatz for the analytic solution of the 
£ — ■ problem in the asymptotic limits of weak and strong coupling. 

The renormalization-group approach to interacting fermions in more than one spatial dimension has been ex- 
tensively applied to the study of instabilities of the Fermi liquid state, and has become a major tool in the study of 
correlated electron systems. Nevertheless, even for weak electron-electron interactions, the picture is far from being 
complete, since electrons in solids also interact with bosonic modes such as phonons. Therefore, the development 
of an RG scheme Q that includes both electron-electron and electron-phonon interactions on an equal footing is an 
^2 . important advance. Experimental evidence show that in many strongly correlated systems, such as high-temperature 
superconductors and organic charge transfer salts, both electron interactions and phonons seem to play an important 
*J i roleiHia. 

The renormalization-group approach to interacting fermions coupled to phonons was presented in Ref. pj. This 
approach takes retardation effects and the presence of multiple energy scales fully into account. For a circular 
1 . • Fermi surface, the RG equations predict the onset of the superconducting instability in agreement withe Eliashberg's 
superconducting theory jsj. A large- N analysis, where N ~ Ep/A is the number of patches in the Fermi surface (Ep 
is the Fermi energy and A is the size of the patch) , shows that in this case Eliashberg theory is asymptotically exact 
q ' and Migdal's theorem [jj emerges as a consequence of the 1/N expansion. Here we present a numerical solution of the 
RG equations at T = 0, showing how the couplings, which are functions of frequencies, flow with the RG procedure. 
In the limits of weak and strong phonon coupling, simple analytical expressions are extracted on the basis of the 
^ | numerical solution. 

\q ■ For completeness, in Sec. [I]the RG equations for the self-energy and interaction couplings are derived. In Sec. ITTlwc 
£NJ ' present the numerical results, and the analytical expressions associated with the asymptotic limits. Sec. IIIII contains 
. the concluding remarks. 
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I. DERIVATION OF THE RG EQUATIONS 

The action that describes electron-electron and electron-phonon interactions can be written as S(ip,(f>) = S e (ip) + 
S P h{4>) + S e - P h{vj, 4>) + S e - e (yj), where (f> are bosonic fields, ip are fermionic (Grassman) fields (we use units such that 
7 H=l = k B ), 

S e = ^kA iuJ ~ e k)lpka , (1) 

i J ojkcr 

is the free electron action, ek is the electron dispersion as a function of momentum k, and 



S P h = / 4>\{i& ~ Wa)4> q , (2) 

is the free phonon action where ui q is the phonon dispersion (cr =|, J, is the electron spin, k = {ui, k} and q = {f2, q}, 
where w, f2 are fermionic and bosonic Matsubara frequencies, respectively, and k, q are the momenta). The electron- 
phonon interaction can be written as: 

Se-ph= / 9{<l)^l +q ai } ka{4>q + 4>-q) , (3) 



where g(q) is the electron-phonon coupling constant. The electron-electron interactions have the form: 

1 3 f 

oil/ u(k i ,kz,k 2 ,k 1 )vjl ia ^ k2IJ i)l 3rT ,vj kia , , (4) 
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where u({ki}) is a general spin-independent electron-electron coupling which depends on the momenta and frequencies 
of the electrons, and k 4 = k\ + k 2 — fc 3 due to momentum conservation. 

Since the bosonic action is quadratic in the boson fields, they can be integrated out exactly, leading to an effective 
electron-electron problem with retarded interactions: 

u(k 4) fe 3 , k 2 , kx) = u(k 4 , k 3 , k 2 , h) - 2g(kx, k 3 )g{k 2 ,k i )D{k 1 - k 3 ) , (5) 

where 

D(q) =wJ{lj 2 +cjI) (6) 

is the phonon propagator. Here we consider the case of a circular Fermi surface and anisotropic Einstein phonons, 
with g(k,k') — go and w q = lje- In what follows, it is also convenient to define the dimensionless electron-phonon 
coupling constant 

A = 2N{Q)^/u E . (7) 

In the Kadanoff- Wilson approach to RG the flow equations are obtained by computing the corrections to the 
couplings of the theory as the energy modes in a energy shell between A and A + dA are integrated out. The 
momenta, frequencies, and fields are then rescaled in such a way that the quadratic terms remain unchanged. This 
second step presents a difficulty in the electron-phonon problem. The momenta of the electrons scale differently in 
the directions parallel and perpendicular to the Fermi surface, while the momenta of the phonons scale isotropically. 
In our RG treatment of the electron-phonon problem yj, we use the quantum field theory version of the RG, with 
no rescaling. In this approach, corrections to the vertices of the model are written in term of running, i. e. cut-off- 
dependent, coupling functions. The RG flow equations for these running couplings are then obtained by imposing 
that the vertices are cut-off independent. 

We start with the two-point vertex, which at one loop is given by 

f°° dO. f d 2 V p im + 

r( 2 )(s( w ,k)) = E^,k)- / ^ u(k,k',k,k') (8) 

7_oo 2?r J (a) (2-kY ill - e k / - 2^(S2,k') 

where k = (w,k) and k' — (O, k'), and X^(w,k) is the electronic self-energy which now flows under the RG, and the 
RG parameter £ = log(A /A), so that d£ = —dA/A. The RG equation for the electron self-energy is obtained by 
imposing the condition of renormalizability of the theory, that is, 

dr^ 

where t = log(A /A) and d£ = -dA/A. 

Notice that the interaction that appears in Eq. (JSJ, u(k,k' ,k,k'), only describes the forward scattering channel 
(&4 = fci = k, ks — k 2 = k') which, in the large N limit, does not flow under RG ,1]. Therefore, it can be substituted 
by its unrenormalized value given by Eq. (|3J). The electron-electron part of the interaction, u(k,k',k,k') = uq 
contributes only to the real part of the self-energy and gives a shift in the chemical potential that can be reabsorded 
into the definition of the chemical potential 2]. The electron-phonon part for w > 0, leads to the following RG 
equation for the imaginary part of the self-energy, E (lu): 

d£^ (w) = - L ~ x ^ n ~ "\n-^ m y + Al ' (10) 

where A^ = Aoe _f (where Ao < Ep is the cut-off in the beginning of the RG flow and it is assumed to be much larger 
than the other energy scales in the problem). There is no dependence of on the direction of the momentum k 
because we are considering the case of isotropic Fermi surface, and the dependence on the magnitude of k is irrelevant 
2J. It is convenient to write S" (uj) = [1 — Z(u))\u), and the solution of Eq. (|10() becomes: 

1 f +co dtt 

z t (w) = i + - —\Lu E D(n-u})nFAn) (11) 

w ./__ 2 



where 



2 r°° dW 

Flip) = ~ / — -o » (12) 
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where we have introduce Wt(u)) = ke,/Zt(iS) is the new renormalized cut-off running scale. The dependence of Zyfiui) 
on W is weak (see below) and therefore we can write Z\y(u>) « Zw>(u>). In this case (|12[) can be solved at once: 



tan 



V M 



(13) 



In what follows we will be only interested in the low frequency behavior, in which case Wi{u>) can be safely replaced 
by Wi(uj = 0), the expression for Zt(u) can be easily evaluated from ijTT|) to give 



Z t (u) = 1 - A 



LJ E 



tan 



i fWe+uj E 



(14) 



Notice that indeed the dependence of Zw on W is weak, as assumed previously. In the static limit, u> — > 0, one 
obtains 



Zi(0) = l + \LO E /{W i +LO E ). 



(15) 



So far we have considered the renormalization of the self-energy. The renormalization of the interaction in the 
Cooper channel (fej = — k% and &2 = — &i) can be obtained in a completely analogous way. Since we are considering 
the case of a circular Fermi surface, we can focus on the s-wave component of the BCS interaction v(u>i,uj3) = 
N(0) J dOi/2-K J d9 3 /2-Ku(—k 3l fc 3 , — ki, fci). The RG equation for this interaction can be shown to be 0: 



—v(ui,UJ2j) 



+ OC 



A, 



7T [ w -e;'h] 2 + a; 



(16) 



where the initial condition for the flow is given by v(<Ji,u>3,l = 0) = uq — \lu e D(lui — 0J3). Equation (|16f) can be 
written in matrix equation as: 



<9U 



= UMU 



where 



Uij(£) = v(uji, ujj,£) 
M l3 (£) = AtSij/ftwi-Sifatf + At] 



Formally, we can rewrite (|17|l as 



d\J- 



d( 



M 



since U^^Utr 1 = -QgU -1 . The solution of HJ reads: 

LT 1 ^) =U- 1 (0) + P(^) 

where 



which can inverted to give: 



P(£)= [ d£'M{£') , 
J 



V(£) = [1 + U(0) -P^T 1 -U(0) 



(17) 



(18) 
(19) 



(20) 



(21) 



(22) 



(23) 



Ea. l|23(l allows the study of instabilities of the Fermi liquid state towards superconducting instabilities. The insta- 
bilities occur when one of the couplings diverges under the RG flow. These instabilities happen at some finite energy 
(or temperature) scale £ c at which one of the eigenvalues of the coupling matrix XJ(£ C ) diverges at £ = £ c . Notice that 
the condition for the instability can be written as: 



det[l + U(0) -P(4)] =0 



(24) 
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Hence, the problem reduces to the calculation of the zeros of a determinant or, equivalently, to the problem of finding 
the zero eigenvalue of the matrix [1 + U(0) • P]: 



[1 + U(0) • P(4)] • f = 
where the fi is the eigenvector of the problem. Eq. (|25(l can be written explicitly as: 



/H = - / — 



f/O 

7T 



dA 



1 



[nZ e (Q)} 2 + A 2 



m 



(lu - n) 2 + u \ 



/(«). 



Here we can, similarly to the expression for the self-energy l|llf) . write 

Zt e (u))(f>t c (w) = - 



dw' 
2 



M 



(lu - UJ 1 ) 2 + 0j\ 



(25) 



(26) 



(27) 



where we have defined <^ c (w) = f£ c (ui)/Z^ c (ui). Equations ifHI) and (|2"T|) determine the energy scale W c = Woe~ ic at 
which the renormalization group equations for the scattering in the Cooper channel diverge as one renormalizes the 
problem from high to low energies. Below this energy scale the Fermi liquid breaks down and superconductivity sets 
in. Thus, we can associate W c with the superconducting gap, A. In fact, we will show that the equation for W c gives 
exactly the same result obtained from Eliashberg's theory for strongly coupled superconductors 



II. SOLUTION OF THE RG FLOW EQUATIONS 

Traditionally the superconducting temperature has been calculated directly from Eliashberg's equations Q . Nev- 
ertheless, the formalism presented in the previous section allows the solution of the problem by solving equations l|14() 
and l|23J) . instead. The advantage of such a procedure is clear since it is not necessary to solve integral equations. 

It is interesting to investigate how the coupling matrix evolves under the RG flow in different regimes since it 
provides an insight on how to solve the problem analytically in some asymptotic limits. The simplest case occurs 
when there are no phonons present (as discussed in Ref|2j) where a full analytical solution is possible. For the case of 
A = 0, Eq. H14f> gives Zi = 1 and i|23|) becomes independent of the external frequency and therefore we must have 
4>(u>) — C is a non-zero constant. In this case (|23[) gives: 

f + °° du' , „ 
1 = -u Q I -z~F tc (u ) 

J-oo * 

= -u \n(W /W c ) (28) 

where we have used (|T3)) . Notice that because W c < Wo the above equation only has solution if mo < 0, that is, if the 
unrenormalized electron-electron interactions are attractive to start with. The solution of the above equation gives: 

A = A e 1/u ° , (29) 

where Wo — Ao since Z = 1. Since, in our case we assume mo > 0, no superconducting instability can exist in the 
absence of phonons. 

Instead of only focusing on the solution of the RG equation at the instability, at which point the zero eigenvalue 
condition l|25|l holds, we have also obtained the full solution for the RG flow by solving and (|23[1 numerically. 
Fig. shows the evolution of U(£) with I. Each panel in Fig. n re P resen ts the NxN matrix U(£) at a given £. The 
frequency has been cut-off so that \oj\ < Ao and has been discretized in N = 200 divisions in this interval. In a typical 
flow in the weak-coupling limit (in the case of Fig. ^ A = 0.3, cue = 10, and Mo = 0.1), U flows from the initial 
condition where v(to, —lu) are the largest couplings to the instability point where the couplings that first diverge are 
the v(u>i, ujj) where |ct>i|, \ujj\ < cue- 

We can clearly see from FigQ] that, as the RG develops, the matrix U, which at the beginning of the RG has 
its largest elements along the diagonal, acquires a cross-like shape, at which point its largest matrix elements occur 
near the origin where the frequencies are very small. The region with largest matrix elements occurs within the dark 
circle of size loe around the origin. Moreover, within this circle the matrix elements are essentially constant showing 
that there is very little frequency dependence if the frequencies are smaller than loe- Hence, loe is the scale that 
separates high from low energy physics in the weak coupling regime. Thus, the full frequency dependence of the 
phonon propagator is irrelevant and we can safely approximate: 



(30) 
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Figure 1: Plots of the NxN matrix U at different RG scales £. Here the number of frequency divisions iV = 200, and the value 
of the parameters used are A = 0.3, Ao = 100, loe = W, Uo = 0.1. Panels correspond to I — 0, 2.5, 3, 5, 6.5, 6.9, 7.1, and 7.19. 



Hence, from Eq. (|14|) we have Zi(u>) « 1 + A which is independent of I. 

Even tough the most divergent elements of the matrix U are the ones at small frequencies, it is important to include 
all the matrix elements since the RG equation couples low and high frequency processes. Again, making use of the 
fact that u> E is the energy scale that separates high and low energy physics in this limit, we can propose the ansatz: 



!>0 



if \oj\ < uje 
if \u>\ > uje 



(31) 



where 0o and are unknowns that to be calculated. Substituting i|30|) . 1311) . and the expression for Zg c into 123|) 
one finds: 



where 



[1 + A + (uo - A) h(u E /A c )} 0o + (1 + A) (4 - h{ujE/Ac))cboo = , 
u h(uj E /A c )(j) Q + (1 + A) [1 + (£ c - h(uj E /A c ))u } ^ = , 



, . . 2 f z tan 1 (x) z 2 
h{z) = - / dx ^ , = — $(-z 2 ) 

7T Jo X 27T 



(32) 
(33) 



and where $(x) is a LerchPhi function, = Efclo + V 2 ) 2 - The solution of ^ 

is given by the determinant 

equation: 



1 + A + (u - A) h{uj E /A c ) (1 + A)(4 - h(u E /A c )) 

u h{uj E /A c ) (1 + A)(l + (4 - h{u) E /A c ))) u a 



0. 



or equivalently, 



Amo^ 2 (w b /A c ) - A(l + u + u Qh(uj E /A c ) + (1 + A) (1 + u £ c ) = , 



(34) 



(35) 



which is a transcendental equation for A c . The last equation can be rewritten in a more interesting form by defining 
A c = x c uj e where x c is the variable of the problem (notice that £ c = ln(A /A c ) = £ E + hi(l/x c ) where £ E = ln(A /w^;)) 
in order to get: 



h 2 {l/x c ) - [1/u* + 1 + ln(l/x c )] h(l/x c ) + + ln(l/.x c )] = . 



(36) 
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where 



l + u ln(A /w E ) 



(37) 



is the renormalized electron-electron interaction at the scale of phonon energy uje- This is the Anderson-Morel 
potential^], which emerges naturally from the solution to the RG flow. 

Note that the renormalized interaction is really weakly dependent on the bare electron-electron interaction since 
as Uq increases it saturates very fast at /i* w l/ln(A /wE) which is interaction independent. In fact, in ordinary 
superconductors /z* varies little from material to material (0.1 < fi* < 0.25) and this can be understood by the 
logarithmic dependence with the cut-off Ao and the phonon frequency loe (if we estimate Ao ~ 10 4 K while cue ~ 10 2 K 
we get fi* « 0.217). 

When x c -C 1 (A c <C lje) we can approximate h(l/x c ) ~ ln(l/x c ) and Eq. l|3l)| becomes 



— + 1 

fj,* 



1 + A 



ln(l/a; c 



1 + A 
A/i* 



which leads to 



W c -loe exp[-(l + A)/(A-/j*)] 
which is the MacMillan expression for the superconducting gap at zero temperature 0, . 



(38) 



(39) 
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Figure 2: Plots of the NxN matrix U at different RG scales £. Here the number of frequency divisions iV = 200, and the value 
of the other parameters are A = 4, Ao = 100, uje = 10, and no = 0.1. Panels correspond to t — 0, 1, 2, 2.5, 3, 3.13, 3.157, 3.172. 
The scale 2W C ~ 40 distinguishes the high and low frequencies close to £ c . 



In the strong-coupling regime (A > 1) the situation is rather different. Fig. shows the numerical results for the 
full solution of the RG flow for parameters in strong coupling. The form of the matrix U(£) is very different from the 
weak-coupling limit and the important scale separating the high and low energy physics is 2W C . We have seen that 
for A < 1 the characteristic energy scale that separates high from low energy physics is wg. When A > 1 we expect 
W c > uje in which case loe is not the characteristic energy scale of the instability. Instead we expect W c to act as the 
characteristic energy scale in separating the low and high energy physics. Thus, it is natural to make the following 
ansatz: 



(40) 
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where 4>o is an unknown. Since, as discussed previously, \i* is bounded from above, the electron-electron interaction 
will play a minimal role in the problem and can be disregarded. 
In this regime (|2*7)l becomes (4>o ^ 0): 



., [* w <> cLu [ + °° IA 
• z 1 1 dA- 



w £ 7 W ffc (w 2 + a 2 )( w 2 + w |) 7T \2w c mi) 



But from O we have Z c (0) ps 1 + \(u)e/W c ) - A(w B /W / c ) 2 , which substituted in (gTJ gives 
the strong coupling result of Eliashberg equation found by Allen and Dynes ^1 • 



(41) 



(42) 



III. DISCUSSION AND CONCLUSION 

We have numerically solved the RG equations for the flow of the BCS coupling matrix for interacting electrons 
coupled to phonons, for the case of circular Fermi surface and Einstein phonons. Fig. |21shows the result for the energy 
scale of divergence in the Cooper channel, for a range of A, and the expected analytical result at weak and strong 
coupling limits. One can clearly see that the approximate analytical solutions, based on the evolution of the coupling 
matrix under the RG, give a very good description of the superconducting gap in the asymptotic limits. Thus, our 
RG scheme can used as an alternative way to solve for the superconducting properties of a Fermi liquid and is in full 
agreement with Eliashberg's theory. 




Figure 3: Plot of W c (~ superconducting gap) in units of uje versus A (black circles) (N = 200, Ao = 100, uje = 10, Mo = 0.1) 
The dashed line is the fit at small A's using the MacMillan formula and the solid line is the fit to yf\ at large A's (Allen-Dynes) 



Unlike mean-held approaches like Eliashberg theory, the RG approach, when solved numerically, can be easily 
extended to any shape of the Fermi surface or anisotropic phonons. It also addresses the question of competition 
between different types of instabilities, as was done extensively for the case without phonons 0,0,0]. Furthermore, 
the approach presented here at zero temperature can be easily extended to hnite temperature. 
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